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Abstract. Observations of clusters of galaxies that gravitationally lens faint background galaxies can probe the 
amount and the equation of state, wx, of the dark energy (quintessence) in the universe. Provided that the mass 
profile and the mass normalization of the cluster are determined, it is possible to constrain the cosmological 
parameters that enter the lensing equations by means of the angular diameter distances, by locating (either by 
observations of giant arcs and magnification bias effect) the critical lines corresponding to known redshift source 
populations of galaxies. This method can help to distinguish between accelerating and decelerating models of the 
universe. Furthermore, since the position of critical lines is affected, especially in low-matter density universes, by 
the properties of quintessence, the observations of a suitable number of lensing clusters at intermediate redshifts 
can determine the equation of state. A very preliminary application of the method to the cluster CL 0024+1654 
seems to support a flat accelerating universe dominated by dark energy. 
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1. Introduction 

Observational cosmology has devoted large efforts in the 
last years to characterize the energy content of the uni- 
verse. Galaxy clustering (Bachall & Fan 1998| ; Carlberg et 



universe (Harun-or-Rashid & Roos |2001| )(»a/ + f^x = 1, 
where fix is the dark energy density parameter of the 
universe) with 0.3-0.4 and a constant equation of 



al. 1998) and large-scale structure (Peacock et al. 2001 



Verde et al. 2001) observations favour models of a uni- 



state -1 < ~ -0.4 (Waga & Miceh |1998| Wang et al. 
[2000 ) at the 68% confidence level or better according to a 
concordance analysis (Wang et al. 2000). A less conserva- 



verse with a subcritical matter energy density ^Im (Turner 
200C). Since, according to balloon-based measurements 
of the anisotropy of the Cosmic Microwave Background 
Radiation (de Bernardis et al. |2000| ; Balbi et al. |2000D , the 
total of energy content of the universe nearly equals the 
critical density (Jaffe et al. p001| Pryke et al. p002[ ), we ex- 
pect that about 2/3 of the critical density is in form of dark 
energy (also called quintessence). Furthermore, evidence 
coming from type la supernovae th at th e universe is accel- 
erating its expansion (Riess et al. 1996 ; Perlmutter et al. 



1999a ) demands a strongly negative pressure for the dark 
energy (wx = Px / Px < —1/3, where px and px are, re- 
spectively, the pressure and energy density of the dark en- 
ergy). These observations, together with other constraints 
coming from the age of the universe, gravitational lens- 
ing statistics and Lya forest, support a geometrically flat 
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five maximum likelihood analysis suggests a smaller range 
for th e equa tion of state, — 1 < wx ^ —0.6 (Perl mutte r 
et al. |l999h| ; Wang et al. |2000| ; Bean & Melchiorri ^002] ). 
According to these results, in what follows, without being 
explicitly stated, we will assume a flat universe. 

After the first proposal of dark energy (the cosmologi- 
cal constant, wx = — 1), many other candidates have been 
suggested. One interesting idea is that the energy density 
is provided by a scalar field rolling down an almost flat 
potential (Caldwell et al. |199^ ; Ratra fc Pee bles |1998[ de 
Ritis et al. ^000| Rubano & Scudellaro [2002] ) . Other pos- 
sibilities are represented by a fluid with a constant equa- 
tion of state, called X-matter (Chiba et al. 1997; Turner & 
White 1997), or by a network of ligh t non-intercom muting 
topological defects (Vilenkin |l984| ; Spergel & Pen |l997t 
{wx = —m/Z where ni is the dimension of the defect: for 
a string, m — 1] for a domain wall, m = 2). Generally, the 
equation of state wx evolves with the redshift, and the 
feasibility of reconstructing its time evolution has been 
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investigated (Cooray & Huterer 1999; Chiba & Nakamura 
Saini et 



200C 



al. 



2000 



Goliath et al. 2001: Huterer & 



Turner |200l| ; Maor et al 2001| ; Nakamura & Chiba 2001 
Wang & Garnavich 200l[ Yam amoto & Futamasc 2001 
Corasaniti & Copeland 2002). Since in flat Fricdmann 



Lemaitre-Robertson- Walker (FLRW) models the distance 
depends on wx only through a triple integral on the red- 
shift (Maor et al. 2001), wx{z) can be determined only 
given a prior knowledge of the matter density of the unl 

Weller 



verse (Gohath et al. 2001 



Gerke 



Albrecht |2001| ; 
& Efstathiou 2002). In what follows, we will consider only 
the case of a constant equation of state. 

Although the listed results are really compelling, it is 
still useful to develop new tools for the determination of 
the cosmological parameters. Many of the discussed meth- 
ods are affected by shortcomings, like poorly controlled 
systematic errors or large numbers of model parameters 
involved in the analysis. An independent constraint can 
improve the statistical significance of the statement about 
the geometry of the universe and can disentangle the de- 
generacy in the space of the cosmological parameters. 

Gravitational lensing systems have been investigated 
as probes of dark energy. Gravitational lensing statistics 
(Wagafe Mi celi |1998[ Cooray & Huterer |1999| ; Wang et 
al. ^000 ; Zhu 200C), effects of large-scale structure growth 
in weak lensing surveys (Benabed & Bernardeau ^001 ) 
and Einstein rings in galaxy-quasar systems (Futamase 
& Yoshida 2001; Yamamoto & Futamase 2001) are very 
promising ways to test quintessence. Here, we propose to 
investigate clusters of galaxies acting as lenses on back- 
ground high redshift galaxies. The feasibility of these sys- 
tems to provide information on th e universe is already 
know n (Paczyhs ki fc Gorski 1981 ; Breimer & Sanders 



1992| ; Fort e t al. |l997| ; Lin k fc Pierce |l99S| ; Lombardi & 



Berlin |l999t Gautret et al. |2000|) . Provided that the mod- 
eling of the lens is constrained, once both arc positions 
and its redshift are measured, it is possible to gain an in- 
sight into second-order cosmological parameters contained 
in angular diameter distances ratios (Chiba & Takahashi 
2002 ; Golse et al. 2002 ). In addition to observations of 
arcs, a statistical approach based on magnification bias 
(Bro adhurst et al. |l995| ; Fort et al. |l997| ; Mayen & Soucail 
200C ) can as well locate the critical lines (locations of max- 
imum amplification) corresponding to background source 
populations. 

In this paper, we will explore the feasibility of clusters 
of galaxies acting as lenses in probing both the amount 
and the equation of state of quintessence in the universe, 
assumed to be flat. In Sect. 2, we outline the method. 
Section 3 is devoted to an application to the cluster of 
galaxies CL 0024-f 1654. In Sect. 4, we discuss some sys- 
tematics affecting the method. Some final considerations 
are presented in Sect. 5. 



2. How critical lines depend on dark energy 

The study of critical lines in a gravitational lensing system 
is a potentially important tool to probe the content of 




Fig. 1. The ratio of distances Dds/Ds as a function of the 
source redshift for a deflector at Zd = 0.3, for different sets 
of cosmological parameters. The thick lines correspond to 
Q,M = 0.3; the thin lines to = 0.5. The full and dashed 
lines correspond to, respectively, wx = —1 and wx = 
-1/3. 
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Fig. 2. Contours of equal Dds/Ds on the {^m, wx) plane 
for Zd — 0.3 and Zs = 1. Each contour is drawn with a step 
of 0.01. The value of the contours increases from the top 
right corner to the bottom left corner. The thin dashed 
lines correspond to lines of constant q^. 



dark energy in the universe and to constrain its equation 
of state, as already shown in the case of galaxy-quasar 
lensing in Futam ase & Yoshida ( ^001 ) and Yamamoto & 
Futamase (2001). This type of cosmological investigations 
requires an accurate modeling of the lens, the observation 
of a critical line and the knowledge of the redshifts of both 



the lens an d the deflected sourc e (Bre imer & Sanders 1992 
Fort et al. |1997|; Link & Pierce |1998[). 
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Fig. 3. The derivative of the ratio of distances Dds/Ds 
with respect to wx for a lens at = 0.3 as a function of 
the source redshift, for different values of the equation of 
state. The full lines correspond to wx = — 1; the dashed 
lines to wx = —2/3; the long-dashed lines to wx — —1/3. 
The thick (thin) lines are for Hm — 0.3(0.5). 
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Fig. 4. The relative variation between the ratio of dis- 
tances Dds/Ds for two cosmological models, (Qm = 
0.3, Wx — —1) and = 0.3, wx = —1/3), in the {zd, Zs) 
plane. Each contour is drawn with a step of 0.01. 



As an example for our quantitative considerations, 
let us consider as deflecting cluster a singular isothermal 
sphere (SIS). The projected density mass E of the SIS is 



1 



2G DdO' 



(1) 



where G is the Newtonian constant of gravitation, a the 
velocity dispersion, Dd the angular diameter distance to 
the deflector and the angular position in the sky. For a 



spherically symmetric lens, the tangential critical line is 
determined by (Schneider et al. 1992) 



UGM{et) D 



ds 



DdDs 



(2) 



where c is the velocity of the light and M{0) is the lens 
mass within the radius 0] Dds and Ds are the angular 
diameter distances to the source from, respectively, the 
lens and the observer. For the SIS, Eq. (0) reduces to 



= 47r 



/cr\2 Dds 
\c) D' 



(3) 



Once 9t and a are known, the ratio of distances Dds/Ds 
can be determined. 

The dependence on the cosmological parameters is 
contained in the angular diameter distance. In a flat 
FLRW universe, the angular diameter distance between 
an observer at Zd and a source at Zs is 



D{zd,Zs) 



1 



Hol + Zs 



(4) 



dz 



. V^Mii + zf + (1 - r!M)(i + z)3(--+i) 



Hq is the today Hubble parameter. At high redshift, the 
pressureless matter density overcomes the dark energy; 
for large Zd and small wx , Dds is nearly insensitive to the 
equation of state. For the expression of the distance in 
inhomogcneous universes, we refer to Sereno et al. ( ^001 
^0021 ). 



Let us go, now, to examine the feasibility of deter- 
mining Wx with observations of strong lensing events in 
clusters of galaxies by the study of the ratio of distances 
Dds/Ds- Once the lens redshift is fixed, Dds/Ds first in- 
creases rapidly with the source redshift and, then, for Zs 
gre ater t han 2.5, is nearly constant (Asada 1997 ; Fort et 
al. 1997), as can be seen in Fig. (^. The change with the 
cosmological parameters can be significant. The ratio in- 
creases with decreasing Qm and with dark energy with 
large negative pressure, i.e. it is maximum in the case of 
the cosmological constant. The variations with r2j\/ and 
Wx are comparable. Changing Qm from 0.3 to 0.5 has 
the same effect of increasing wx from —1 to —1/3, so 
that Dds/Ds is nearly indistinguishable in a universe with 
= 0.3 filled in with string networks and in a model 
with = 0.5 and cosmological constant. 

To quantify the dependence of Dds/Ds with the cosmo- 
logical parameters, we consider fixed redshifts for the lens 
and the source, see Fig. (|^). The ratio is quite sensitive to 
Hm- The variations due to changes in f^M for wx = const. 
are greater than in the case of the constant deceleration 
parameter qo = (l-|-3wx(l — f^M))/2. For Zd = 0.3, Zs = 1, 
and ranging from to 0.6, when ~ Q the variation 
is ~ 4%; when wx — —1, the variation is ~ 15%. The 
dependence on the cosmological parameters is maximum 
for high negative values of go, i.e. the region today pre- 
ferred by observations. For some particular pairs {zd,Zs), 
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i.e. for low lens redshifts and sources very near to the de- 
flector, the ratio is nearly constant on lines of constant 
deceleration parameters; these properties suggest that the 
method of the critical line can help to distinguish between 
accelerating and decelerating universes. The dependence 
of Dds/ Ds on the equation of state increases for low mat- 
ter density universes and the sensitivity nearly doubles for 
small changes in VLm'- for Zd = 0.3 and Zs — 1, the relative 
variation from wx — —1 to wx — ^1/3 is 1.9% (3.4%) 
when r^A/ = 0.5 (0.3). The sensitivity is maximum for in- 
termediate Wx] for large negative pressure {wx ~ -0.9), 
the ratio is nearly independent of variations of the equa- 
tion of state. 

In Fig. (^), the derivative of the ratio Dds/Dg with 
respect to wx is plotted as a function of the redshift of 
the source once the redshift of the deflector is fixed. The 
derivative is negative for a large range of redshifts of both 
source and deflector. Transitions from negative to positive 
values occur for very negative wx- The source redshift 
where the derivative cancels out decreases with increasing 
r^M and Zd'- for Zd = 0.3 (0.6), SIa/ = 0.3 and wx = —1, 
the derivative is null at Zs ~ 7.0 (2.1). The sign of the 
derivative determines, when the equation of state changes, 
as the angular position of the critical lines moves: when the 
derivative is negative (positive), as the equation of state 
increases (i.e. as wx moves from —1 to 0), the angular 
radius in the sky of the critical line, for fixed source and 
deflector redshifts, decreases (increases). 

The modulus of the derivative is an estimate of the 
dependence of the ratio on wx ■ Independently of the value 
of Zd^ and wx^ the dependence on wx first increases 
and takes its maximum at an intermediate source redshift, 
and then decreases quite slowly. For dark energy in the 
form of a cosmological constant {wx ~ ^1); = 0.3 and 
l^Af = 0.3, the maximum is at Zs ^ 0.75. For increasing 
Wxi the maximum moves to higher redshifts: for domain 
walls {wx — —2/3), the maximum is at Zg ^ 1.02. From 
Fig. (^, we see that for a large range of wx and ^Im 
the maximum is at Zs ^ 2. This trend of the derivative 
is connected to the properties of the ratio Dds/Ds, that 
flattens at higher source redshifts. 

Now, we want to search for the optimal lens and source 
configuration in order to discriminate among quintessence 
models. For illustration, we choose two universes with the 
same content of matter {^m = 0.3) but different wx', we 
consider a cosmological constant {wx = —1) and string 
networks {wx — —1/3). In Fig. (P, we scan the {zd,Zs) 
plane plotting the relative variation between the two pairs 
of cosmological parameters. For a given lens redshift, the 
best Zs is very close to the deflector, i.e. a couple of red- 
shifts corresponding to the rising part of the ratio Dds/Ds] 
the sensitivity decreases for larger and larger source red- 
shifts. So, the configurations with high sensitivity to the 
quintessence are those with very low cross section for 
strong lensing events. On the other hand, given a back- 
ground population at z^ ^ 1, the optimal lens is a quite 
high redshift cluster at Zd ^ 0.7; however, the depen- 



dence on the quintessence is nearly constant for lenses at 
Zd ^ 0.6 

In order to estimate the accuracy of the determination 
of the equation of state, the variation induced on Dds/Ds 
by Wx must be compared to the error within which the 
parameters of the lens are known. For the SIS, the error 
in the estimate of the ratio of distances is 



A 



Dds 
Ds 



Aa 
a 



D. 



ds 



Ds 



(5) 



where Act and A9t are the errors, respectively, on the ve- 
locity dispersion and 0t. Usually, the largest uncertainty 
in the modeling of a lens comes from the error in the mea- 
surement of the velocity dispersion. Catalogues of galaxy 
velocities in lensing clusters are of the order of 50, so that 
the uncertainty on ct is ~ 15%. A6t conies from the accu- 
racy of the location of the arc and its radial thickness and 
from the uncertainty on the geometrical properties of the 
lens, i.e. the accuracy of the location of the center, typi- 
cally chosen to coincide with the brightest cluster galaxy, 
and the ellipticity of the mass distribution. For tangen- 
tial arcs at 9t ~ 20" (Williams et al. 1999), an error as 
large as ^ 1" can contribute a 5% error. The error on dt is 
generally negligible with respect to the error in the mass 
normalization and will not be considered in the rest of this 
section. The variation on Dds/Ds connected to changes in 
the equation of state can be expressed as 



d 



dwx 



Dd^ 
Ds 



Awx, 



and so, comparing Eq. (|^) and Eq. 
have a statistical error of 



Awx 



(6) 



for N clusters we 
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Act 




V 




CT 
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Act 




d 




CT 


( 


dw 




(7) 



where the average is on the redshifts of the critical lines. 
The error in the determination of wx increases with Act 
and decreases with the derivative. Since the error induced 
by the velocity dispersion is proportional to the ratio of 
distances Dds/Dg, see Eq. (|]), and the variation induced 
by Wx is proportional to the derivative, see Eq. (^), the 
uncertainty in the estimate of wx is inversely proportional 
to the logarithmic derivative of Dds/Ds, i.e. to the relative 
variation of Dds/Dg. The properties of the logarithmic 
derivative with respect to the cosmological parameters fi^f 
and Wx are the same of the ordinary derivative; the main 
difference is the disappearance of the minimum. As we 
have seen before, the uncertainty in the equation of state, 
given a deflector redshift, increases with Zs and decreases 
for quintessence with wx far away from —1. 

The case of the cosmological constant is the more prob- 
lematic one since the derivative can cancel out (when 
VLm = 0.3 and Zs = 1.5, the derivative is null at Zd — 1.13). 
However, clusters at intermediate redshift {zd ^ 0.4) are 
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quite stable with respect to the error in the equation of 
state. 

As we shall see in the next section, it is possible to ob- 
tain information from a single cluster of galaxies on more 
than one critical line. So, using in Eq. (j^) the number 
N of clusters, the lower limit on Awx is overestimated. 
Given a typical error of '--^ 15% on cr, we can use Eq. (Q) to 
estimate the number of deflectors necessary for estimat- 
ing wx within a given uncertainty. For mean redshifts of 
(zd) — 0.4 and (zg) — 1.2, an uncertainty of Awx — 0.25 
needs ^ 75 (~ 120) lensing clusters in a universe with 
riM = 0.3 and wx = —1/3 (—0.5). N increases with dark 
energy with large negative pressure and large values of 
fijv/- As discussed, the method is unable to constrain the 
equation of state in the extreme case of a cosmological 
constant, when Awx — 0.25 needs ~ 800 clusters and 
Awx — 0.5 needs ~ 200 clusters. In general, to distin- 
guish dark energy with an intermediate value of wx from 
a cosmological constant at 95% confidence level, in a low 
matter density universe, we need 100-200 strong lensing 
events. These simple estimates are in agreement with the 
results in Yamamoto & Futamase ( ^001 ). 

Together with spectroscopic analyses. X-ray observa- 
tions of a lensing cluster can help to estimate the absolute 
mass of the deflector. The projected X-ray cluster mass, 
under the hypotheses of isothermal and hydrostatic equi- 
librium, is proportional to the cluster gas temperature, 
Tx, and Dd (Wu p000| ): X -ray data alone cannot deter- 
mine the mass without a prior knowledge of cosmological 
parameters. However, it has been shown that the relation 
between a and Tx is not affected by cosmic evolution and 
is co nsisten t with the isothermal scenario, a oc T^^ (Wu 



et al. |1998[ ) . Once calibrated this relation. X-ray observa- 
tions obtained with the new generation of telescopes can 
considerably enlarge the data sample of lensing clusters 
with known mass and help to disentangle the effect of 
cosmology and mass normalization of the deflector. 

3. CL 0024+1654 

Now, let us consider the application of the method out- 
lined in Sect. 2 to a well studied cluster of galaxies, 
CL 0024-1-1654, in order to test the feasibility of what 
we are proposing, and how good the results can be. 

CL 0024-1-1654 is one of the best investigated lenses 
in the universe. It is an optically rich cluster of galaxies, 
with a relaxed structure without a single central dominant 
cluster galaxy, at z = 0.395 and with a velocity dispersion 
of (T = 1050±75 km s "^ (D ressier et al. |1999| ; Czoske et 
al. ^OOll Czoske et al. [20021 ). This is the formal velocity 
dispersion estimated with the assumptions of virial equi- 
librium and random galaxy velocities, so that the reported 
error is a purely statistical one. We will consider the ef- 
fect of some possible systematica in the next section. This 
value o f a is consistent with lensing obs ervat ions (Shapiro 
& Iliev^OOOl ). X-ray data (Soucail et al. ^OOOj Bohringer et 
al. 200C ) also support a regular morphology with no signif- 

+4.9 
-2.1 



is compatible with the observed velocity dispersion. A sin- 
gle background galaxy behind CL 0024-1-1654, at spectro- 
scopic redshift z — 1.675 (Broadhurst et al. 2000), is im- 
aged in a well known multiple arc at dt = 30.5" (Kassiola 
et al. |1992| WaUington et al. |1995|, Small et al. p996| Tyson 



et al. 1998). Images are characterized by a bright elon- 
gated knot, surrounded by a low surface brightness halo. 
The knot comprises two peaks, with separations ranging 
from 0.5" to 1.1", roughly consistent with the relative 



lengths of the various arc components (Small et al. 1996). 
Given this peculiar morphology, we assume an indetermi- 
nation on the critical radius A9t ~ 0.7". We do not take 
into account the error on the position of the centre; in the 
analyses considered here, it is determined as a free param- 
eter in the lensing reconstruction. Based on deep images 
with the Hubble Space Telescope, Tyson et al. (199S) per- 
formed a multi-parameter fit, including a number of small 
deflecting "mascons", to the mass profile. Each mascon 
was parameterized with a power-law model (Schneider et 
al. |1992D, 



^{9) cx 



1 + 7 f 



1 



2-7 



(8) 



where 9c is the core radius and 7 is the slope. Remarkably, 
they found that more than 98% of the cluster matter 
is well represented by a single power-law model centred 
near the brightest cluster galaxies with 9c = 10.0±0.9" 
and 7 = 0.57±0.02, slightly shallower than an isothermal 
sphere (75/5 — 0.5). To disentangle the effect of cosmology 
and absolute mass, we have to fix the central den sity S p 
independently of lensing data. It is (Bonnet et al. 1994 ), 



So 



T 2 

iaa ^(i+a)/2 1 



3ttG9c 



Dd' 



(9) 



with a — 2(1 — 7) ^^^d — {l + u^)^'^du. For a power- 
law model, the angular position 9t of the tangential critical 
line is related to the angular diameter distances and the 
parameters of the lens by 



(10) 



Substituting in Eq. ( p^ for Sq and using the fit parame- 
ters, we get an estimate for the ratio Dds/Dg- It is 



D, 



Zd = 0.395, = 1.675) 



76+018 

^- '"-0.12- 



(11) 



leant substructures. The measured value of Tx = 5.7 



The main term in the error budget comes from the inde- 
termination in the velocity dispersion which contributes 
~ 75% of the total error. In Fig. (||), we show the de- 
pendence of Dds/ Dg on the cosmological parameters for a 
lens-source configuration as in CL 0024-1-1654; the values 
of (r^A/, Wx) compatible with the estimate in Eq. ([ll| ) are 
also plotted. Low matter density universes {^M ~ 0.55), 
which are accelerating their expansion, are favoured. We 
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find —1 < wx ~ —0.2, with the lower values of Hm cor- 
responding to the higher values of the equation of state. 

The method of the depletion curves, i.e. the variation 
along the radial direction in the surface density of back- 
ground galaxies around a massive cluster of galaxies, has 
be en em ployed to further study CL 0024-1-165 4 (Fo rt et 
al. |l997|; van Kampen |1998| ; Rognvaldsson et al. |200l| ; Dye 
et al. ^002 ). Observations of the magnification bias have 
been obtained in the B- and /-band (Fort et al. 1997 ) and 
in the U - and i?-band (Rognvaldsson et al. 2001; Dye et 
al. 2002 ). Extrapolating Hubble Space Telescope data to 
their detection limit. Dye et al. (2002) obtained, for the 
background i?-galaxies, a mean redshift of (zg) — 1.2±0.3. 
From a fit to the SIS profile of the depletion curve in the 
i?-band, the location of the critical curve comes out at 
15±10" (Dye et al. |2002| ). Using these estimates in Eq. (|), 
we can obtain a second constraint on the ratio Dds/Ds] 
unfortunately, as can be seen in Fig. (^ , the uncertainties 
completely hide the second order effect of the cosmological 
parameters on the ratio of distances Dds/Ds- 

A more interesting result can be obtained from the 



/-band. As discussed in Fort et al. (1997), the angular ra- 
dius where the depletion curve starts to increase locates 
the last critical line, that is the critical line corresponding 
to the farther source population. The last critical line at 
^ 60" in the /-band corresponds to background galaxies 
at redshift 2.5 < z < 6.5; however, about 20% of the very 
faint / selected g alaxie s should be above z = 4. As noted 
in van Kampen (199S), given the very low density of the 
background /-galaxies, an appropriate radial binning to 
study the radial profile of the magnification bias is 30". 
So, we will consider an error of 15". This estimate of the 
location of the last critical line is independent of the as- 
sumed mass profile, and can be used in Eq. ( |l0| ) to obtain 
a new constraint on Dds/Ds, see Fig. (^. Since Dds/Dg is 
nearly flat for Zg ~ 2.5, the value of the ratio of distances 
is quite insensitive to the value of Zs corresponding to the 
last critical line. 

Some interesting considerations are obtained from the 
variation of the ratio Dds/Ds with the redshift of the 
source. Figure (^ shows Dds/Ds for a lens at z = 0.395 
and for various cosmological models. Models without dark 
energy are rejected, with no regard to the value of the 
pressureless matter density: both open (iljvf < 1) and flat 
(the Einstein-de Sitter model, with Qm = 1) dark matter 
models are very poorly consistent with the experimental 
points. On the other hand, flat universes with quintessence 
are in agreement with the data. In particular, the data 
from the /-band analysis are marginally compatible (at 
the 68% confidence level) with a flat de Sitter universe 
{^x = 1 and Wx = — 1). Given the large uncertainties, 
we cannot draw definitive conclusions on this multi-band 
analysis. However, even if the data from the /^-galaxies 
have not information on the cosmology, the data from the 
multiple arc and the last critical line in the /-band prefer 
accelerating universes with subcritical matter density. 
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Fig. 5. Contours of equal Dds/Ds on the {flM, wx) plane 
for CL 0024-1-1654 {zd = 0.395) and its multiple arc {zs = 
1.675). Each contour is drawn with a step of 0.01. The 
value of the contours increases form 0.61 in the top right 
corner to 0.76 in the lower left corner. The thick lines 
correspond to the data in Eq. (pl|). The thick full line 
corresponds to the best parameters; the dashed one to the 
lower limit. The thin dashed line separates accelerating 
universes (below) from decelerating ones (above). 




Fig. 6. Dds/Ds as a function of the source redshift for 
CL 0024-1-1654 {zd = 0.395) for different cosmologi- 
cal models. The thick lines are for flat models with 
quintessence. The full thick line for D,x = 1 and wx = — 1 
(de Sitter universe); the dashing thick line is for — 0.3 
and Wx = ^1- The thin lines are for universes with pres- 
sureless matter alone. The full line is for an Einstein-de 
Sitter universe (J^m — 1); the thin dashed line is for an 
open universe {^m — 0. 1). " R-6and" indicates the data 
derived from Dye et al. ( ^002] ) (the error in a is not con- 
sidered); "ARC" is the data in Eq. ([ll]); ^^l-band" indicates 
the data from the depletion curve in Fort et al. ( |l997D . 
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4. Systematics 

In the previous section, we performed a statistical analysis 
based on the data found in the literature. We want now 
to address some systematics that can affect our results. A 
very accurate knowledge of the mass distribution of the 
lens is required to put meaningful constraints on cosmo- 
logical parameters. One of the more important source of 
indetcrmination comes from the modeling of the mass pro- 
file of the lens (Chiba & Takahashi [20021 ). In Eq. (|), we 
have considered only the error coming from a not very ac- 
curate mass normalization but, in general, we have also 
to face the indetcrmination on the cluster mass profile. 
As a general feature, the three-dimensional mass density 
of a clump, p, is proportional to a typical length scale, 
Tg, so that, with respect to the angular diameter distance, 
p cx -D^^. The mass enclosed within an angular radius 9 
comes out 

M{e) cx Daa^Ve, 

where Dd contains the cosmological dependence and 
stands for an overall normalization. 7^ is a factor account- 
ing for the deviations of the cluster mass profile from the 
SIS; 7^ is a function of and of some parameters, such 
as a core radius. Substituting in Eq. (^, we get, for a 
spherically symmetric lens, 

Dels 9t 

To consider the uncertainty in the profile, we have to add 
in quadrature to the right hand side of Eq. (|^) an addi- 
tional relative error of AV/V. Usually, the main contri- 
bution to the error budget comes from the mass normal- 
ization but, in some extreme cases, the indetcrmination 
in the mass profile can be of the same order. AV is max- 
imum when calculated between a lens with a soft core 
and a halo with a singular density steeply rising towards 
the center, as predicted by numerical simulations in the 
standard cold dark matter framework of structure forma- 
tion and approximated by a Navarro- Frenk- White profile 
(NEW, Navarro et al. ( |l995i [l996[ )). 

A NEW model can match the mass distribution of 
CL 0024-^1654 (Broadhurst et al. |2000[). The required 



mass inside the arc's radius for such a model, reproduc- 
ing the projected mass distribution outside the core ra- 
dius, is 40% higher than the prediction of a power-law 
model (Tyson et al. 1998). Without an independent in- 



formation, AV would be a really large error. Eortunately, 
the NEW profile is discarded since it implies a velocity 
dispersion much higher than the measured value (Shapiro 
& Iliev 2000). Once we can discard models with singular 
central density, the indetcrmination in the model can be 
accounted for by the errors within which we constrain the 
parameters of the mass profile (A7 and A9c in the power- 
law model used in Sect. 3). In general, an uncertainty in 
the cluster mass profile can significantly weaken the re- 
sults on the cosmological parameters; but, in the case of 
CL 0024-1-1654, the degeneracy in the fit is not the main 
error. 



Together with the overall mass profile, sub-structures 
must be considered. In the case of a lens with a rather reg- 
ular morphology, even if a "not correct" potential shape 
is used in the reconstruction or the contribution of small 
sub-structures is neglected, the cosmological parameters 
are still retrieved, although with larger errors (Golse et 
al. 2002). On the contrary, neglecting a sub-structure as 



large as 20% of the total mass in a bi-modal cluster com- 
pletely hides the effect of cosmology (Golse et al. 2002). 
Adding the contribution of individual galaxy masses is 
also useful to tighten the confidence levels and can be- 
come critical in some extreme cases, as a galaxy strongly 
pertu rbing the location of multiple-images (Golse et al. 
^002| ). Deep imaging of CL 0024+1654 has made it pos- 
sible to construct a high-resolution map of the projected 
mass distribution of the cluster and to take into account 



the effect of perturbin g gala xies. Kassiola et al. ( 1992 ) 
and Wallington et al. (1995) considered the perturbing 
potentials of two galaxies near the middle segment of the 
arc. Tyson et al. (1998) assigned one or more mascons to 
each of the 118 cluster galaxies and 25 free mascons for 
the remaining cluster mass. However, all these studies in 
literature agree on a overall representation as the one in 
Eq. (1). 

Some features in the 3-D space, as a possible merger 
scenario, can invalidate our estimation of the cosmological 
parameters. A recent analysis of the distribution of the 



galaxies in the redshift space (Czoske et al. 2001, 2002) 
suggests a fairly complicated structure. A group of galax- 
ies lying just in front of the main cluster could be the 
result of a high speed collision of two smaller clusters 
with a merger axis v ery nearly parallel to the line of sight 
(Czoske et al. 2002 ). In particular, a bulk velocity com- 
ponent present in the central velocity distribution would 
over-estimate the mass obtained from the formal central 
velocity dispersion. Eurthermore, galaxies at large pro- 
jected distance from the centre are also affected by the 
collis ion and cannot be used to derive a (Czoske et al. 
^002 ). The consequences on the dark energy constraint 
are quite dramatic since such a scenario could entail a 
systematic error on the estimation of the velocity disper- 
sion of the same order of the statistical one. As can be 
seen from Eig. (^, this additional error would completely 
hide any dependencies on the cosmological parameters. 

Even if, together with strong lensing observations, 
both a weak lensing analysis out to 10 arcmin (Bonnet 
et al. 1994) and X-ray observations (Soucail et al. ^000 



Bohringer et al. 2000) favour a regular morphology, the 
points just discussed suggest caution in the interpretation 
of the results obtained in the previous section. 

5. Some final considerations 

We have explored the feasibility of reconstructing the 
properties of the dark energy in the universe by using 
strong lensing systems in which a cluster of galaxies acts 
as deflector. With respect to other lensing systems, for the 
one discussed in this paper it is possible to determine the 
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position of the critical lines in two independent ways: with 
giant arcs or with the radial shape of the depletion curves. 
This circumstance allows to study the ratio of angular di- 
ameter distances that characterizes the angular position of 
the critical lines over a large range of source redshifts, just 
for a single lensing cluster. Provided that the properties 
of the background populations are well constrained, it is 
possible in principle to use multi-band depletion measure- 
ments to obtain several independent estimates of the ratio 
of distances, each one probing a different source redshift. 

For a flat universe, the sensitivity of the angular posi- 
tions of the critical lines on quintessence becomes higher in 
low-density pressureless matter universes and for dark en- 
ergy with intermediate equation of state. While the anal- 
ysis of only a few lensing clusters suffices to distinguish 
between accelerating and decelerating models of universe 
(also without a prior knowledge of I^m), a considerably 
larger sample (TV ^ 200) and an accurate estimate of flM 
are needed to constrain the equation of state within an 
uncertainty of Awx ~ 0.25 and discriminate, at the 95% 
confidence limit, between a cosmological constant and an 
evolving quintessence. 

In our opinion, a first application of the method to 
the cluster CL 0024+1654 has given interesting results. A 
combined analysis of both the multiple arc and the deple- 
tion curves disfavours models of a universe without dark 
energy. On the other hand, flat accelerating universes are 
in agreement with the data. These preliminary estimates 
agree with the currently favoured constraints from other 
independent measurements. However, some features in the 
redshift space of CL 0024-1-1654, as a possible merger sce- 
nario, could invalidate our results. Indeed, a very accu- 
rate knowledge of the absolute mass distribution of the 
defiector and a correct understanding of the pattern of 
sub-structures are necessary to obtain secure constraints 
on the cosmological parameters. 

The method we have discussed is quite general and can 
be applied to several strong lensing systems. For example, 
a single galaxy, whose stellar velocity dispersion can be 
accurately measured, can multiply image a background 
quasar. Clusters of galaxies need an accurate modeling of 
the pattern of substructures and present a quite problem- 
atic measurement of a but allow one to study the ratio 
Dds/Ds at different source redshifts. Furthermore, a mul- 
tiple image system of galaxies with known redshift makes 
possible an absolute calibration of the total mass of the 
cluster (Golse et al. 2002). We have shown how analyses 
of magnification bias in multi-band photometry can be 
combined with observations of giant arcs to obtain some 
insight on cosmological parameters. 
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